Stationary states of a rotating Bose-Einstein condensate: routes to vortex nucleation 
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Using a focused laser beam we stir a 87 Rb Bose-Einstein 
condensate confined in a magnetic trap. We observe that 
the steady states of the condensate correspond to an elliptic 
cloud, stationary in the rotating frame. These steady states 
depend nonlinearly on the stirring parameters (amplitude and 
frequency), and various solutions can be reached experimen- 
tally depending on the path followed in this parameter space. 
These states can be dynamically unstable and we observe that 
such instabilities lead to vortex nucleation in the condensate. 
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Superfluidity, originally discovered and studied in the 
context of superconductors and later in the system of 
superfluid liquid Helium, is a hallmark property of inter- 
acting quantum fluids and encompasses a whole class of 
fundamental phenomena 0J^]. With the achievement of 
Bose-Einstein condensation in atomic gases [p| , it became 
possible to study these phenomena in an extremely dilute 
quantum fluid, thus helping to bridge the gap between 
theoretical studies, only tractable in dilute systems, and 
experiments. One striking consequence of superfluidity 
is the response of a quantum fluid to a rotating pertur- 
bation. In contrast to a normal fluid, which at thermal 
equilibrium will rotate like a solid body with the per- 
turbation, the thermodynamically stable state of a su- 
perfluid involves no circulation, unless the frequency of 
the perturbation is larger than some critical frequency, 
analogous to the critical velocity B. Moreover, when 
the superfluid does circulate, it can only do so by form- 
ing vortices in which the condensate density vanishes and 
for which the velocity field flow evaluated around a closed 
contour is quantized. 

In this Letter we present the study of the response of a 
Bose-Einstein condensate (BEC) confined in a magnetic 
trap to a rotating perturbation created by a stirring laser 
beam. We observe that for a given perturbation ampli- 
tude and frequency the steady state of the condensate 
corresponds to an elliptic cloud, stationary in the rotat- 
ing frame as predicted in Q . Depending on the path fol- 
lowed in the parameter space (amplitude and frequency) 
of the rotating perturbation, we show that two steady 
states can exist, corresponding to different ellipticities of 
the cloud. We also observe that these states possess an in- 
trinsic dynamical instability Mpl ; moreover, we find that 



this instability can lead to the transformation of the ellip- 
tic state into a state of one or more vortices. The fact that 
vortex nucleation seems to occur only via a dynamic in- 
stability (apart from the phase printing method explored 
in explains why the frequency range over which vor- 
tices are generated is notably smaller than that expected 
from thermodynamics ( |Q and references within). 

A dilute, interacting Bose gas with a large number of 
atoms is well described at low temperature by the hydro- 
dynamic equations for a superfluid These equations 
were studied in the case of a dilute BEC in a rotating 
harmonic trap characterized by the trap frequencies ujx , 
toy , u z [Q (see also || ) . The axes X, Y rotate at the 
frequency £1 around the z axis. There exist solutions for 
which the BEC wavefunction is given by 
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The condensate density for these solutions is a paraboloid 
(Thomas-Fcrmi profile) 
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where g = ATrh 2 a/m is the collisional coupling constant 
set by the s-wave scattering length (a = 5.5 nm for 87 Rb 
in the \F — 2, nip = 2) ground state), m is the atomic 
mass, ii e is the chemical potential in the rotating frame 
(the density is understood to be vanishing for p < 0). 
The XY ellipticity of the paraboloid is related to the 
parameter a appearing in the phase of Eq. (nj) 
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The value for a is determined by the parameters of the 
rotating perturbation according to 
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where a — a/co, = Q/co, 
trap deformation e = (lo\ - 



= \J {uj 2 x + uj Y )/2, and the 
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to three possible values for a corresponding to stationary 
states of the condensate for a given pair (ri,e). These 
solutions are plotted as lines in Figs. 1 and 2. When 
fi varies with a fixed e (Fig. 1), or when e varies with 
a fixed Cl (Fig. 2), the possible values of a are located 
on two branches. Branch I is a monotonic function of f2 
or e, while branch II exists only for some range of the 
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parameters and exhibits a "back bending" , opening the 
possibility for hysteretic behavior. 

The experimental study of these states relies on the 
realization of a rotating harmonic potential and a mea- 
surement of the density profile of the condensate cloud. 
The atoms are confined in an Ioffc-Pritchard magnetic 
trap which provides a static, axisymmctric harmonic po- 
tential of the form U(r) = muj 2 {x 2 + y 2 )/2 + mu) 2 z z 2 /2. 
The condensate in this potential is cigar-shaped with a 
length to diameter aspect ratio of A = uJt/u z which was 
varied between 10 and 25 {uo z /2tt = 11.8 ± 0.1 Hz). The 
atomic cloud is stirred by a focused 500 fiW laser beam 
of wavelength 852 nm and waist wq = 20/Ltm, whose po- 
sition is controlled using acousto-optic deflectors (?J (see 
also jnj). The beam creates an optical-dipole potential 
which can be approximated by mui 2 (exX 2 + eyY 2 )/2. 
The XY axes of the optical potential are rotated at a 
frequency fl producing a rotating harmonic trap charac- 
terized by the three trap frequencies uj\ y — u t + 
and u> z . The value of the transverse trap frequencies is 
measured directly by determining the frequency domain 
(wy < n < ll>x) in which the center of mass motion of the 
cloud becomes dynamically unstable. This measurement 
provides a determination of e and Q 
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To measure the state of the condensate, characterized 
by a, we perform a field-free expansion of the condensate 
for a duration of 25 ms, followed by resonant absorption 
imaging technique along the stirring axis |ll| . Using the 



formalism presented in |12|, we have checked that the 
value of a for the rotating states described by Eq. (|l|) 
changes by less than 10% during the expansion for the 
parameter range studied. 

The experimental procedure begins with the prepara- 
tion of a condensate in the pure magnetic potential by a 
radio-frequency (rf ) forced-evaporation ramp lasting 25 s. 
Starting with 10 s atoms precooled to 10 /iK in an optical 
molasses, the atomic cloud reaches the critical temper- 
ature T c ~ 500 nK with an atom number of ~ 2.5 10 6 . 
The evaporation is continued below T c to a temperature 
of or below 100 nK at which point 3 (±0.7) 10 5 atoms are 
left in the condensate. The rf frequency is then raised to 
a value 20 kHz above t'™ 1 ", the rf frequency which corre- 
sponds to the bottom of the magnetic potential. This rf 
drive is kept present in order to hold the temperature ap- 
proximately constant. At this point, the stirring laser is 
switched on and the condensate is allowed to evolve in the 
combined magnetic and optical potential for a controlled 
duration after which we perform the field free expansion 
and optical detection. 

Since the condensate is first created in a static har- 
monic potential and then the rotating perturbation is in- 
troduced, the final state of the condensate depends on its 
evolution in a time-dependent rotating potential charac- 
terized by the two parameters e(t) and f2(i). If the time 
dependencies are slow enough, the condensate at every 
instant is in a stationary state corresponding to the in- 



stantaneous values of e and f2. In this case, the initial 
state of the condensate is a = 0, and as e and Ct evolve 
so does a evolve along a path defined by Eq. (0). 
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FIG. 1. Steady state value of a as a function of stirring 
frequency Q. The results of the ascending ramp are shown by 
a filled dot (•) and of the descending ramp by a hollow dot 
(o). The branches of Eq. (||) for e = are shown as dashed 
lines. The best agreement between theory and experiment 
was achieved for all data in Figs. 1 and 2 with e — 0.022 and 
a scaling of lj/2it — 200 Hz, while the measured values are 
e = 0.025 ± 0.005 and lu/2-k = 195 ± 1 Hz. This 2% dis- 
crepancy in the frequency may be due to a deviation from 
the Thomas Fermi limit. The angular momentum of the con- 
densate measured after a 200 ms relaxation time reveals the 
presence of vortices for the ascending ramp above Cl = 0.75. 

The first study that we performed involved switch- 
ing on the optical potential with a fixed value of e = 
0.025 ± 0.005 and an initial frequency of Cl = 0. The 
frequency f2 was increased at a constant rate to a fi- 
nal value in the range (0.5,1.5) The rotation fre- 
quency was then held constant at this final value dur- 
ing 5 rotation periods, and finally the state of the con- 
densate was measured. For each final frequency studied, 
two different measurements were performed: a determi- 
nation of a made according to Eq. (|J) and a measurement 
of the angular momentum of the condensate associated 
with the presence of a vortex. The latter measurement, 
which relies on the excitation of a quadrupolar oscillation 
p4] , p^[ , is performed 200 ms after the stirring anisotropy 
is switched off. This delay is (i) short compared to the 
vortex lifetime [Q and (ii) long compared to the mea- 
sured relaxation time (~ 25 ms) of the ellipticity of the 
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rotating condensate. Hence a nonzero value of L z in this 
measurement is the signature of vortex nucleation. Each 
measurement relies on a destructive resonant-absorption 
imaging, and is performed independently; the results are 
shown in Fig. 1 . We repeated the above study with a de- 
scending frequency ramp starting above 1 and finishing 
in the same range as the ascending ramp, and the results 
are also shown in Fig. 1. 




FIG. 2. Measurement of the steady state value of a ob- 
tained from a linear ramp of the stirring anisotropy for two 
different stirring frequencies, (a) When fl = 0.70 the conden- 
sate follows branch I and no vortices are nucleated. Branch I 
of Eq. (Q) is plotted for lo/2tt — 200 Hz (solid line) and also 
for lj/2tt = 197 Hz and lj/2h — 203 Hz (dotted and dashed 
lines), (b) When Q = 0.75 branch II is followed to the point of 
back bending(e c = 0.023). The solid line represents the solu- 
tions of Eq. (ffl for u}/2tt = 200 Hz. (c) For ramps which pass 
this point the angular momentum measured after a 200 ms 
relaxation time reveals that vortices are nucleated. 

In both cases, the state of the condensate starts on and 
follows one of the branches of Eq. (||). For the descend- 
ing ramp, the lower part of branch II is followed until 
the back bending frequency is reached below which this 
branch ceases to be a solution of Eq. 0. At this point, the 



condensate is observed to switch to branch I and to follow 
it down to SI = 0. The fact that no vortices are nucleated 
along this path seems to confirm the prediction that the 
lower part of branch II is stable Q|. For the ascending 
ramp the condensate is observed to start on and to fol- 
low branch I until a frequency of Sl c ~ 0.75 is reached. 
Beyond this point the experimentally determined value 
of a is no longer a solution of Eq. (^), and the presence of 
vortices is detected in the angular momentum measure- 
ment of the condensate. The value of Sl c is in good agree- 
ment with the frequency at which branch I is expected to 
become dynamically unstable 0]. This experiment was 
repeated for different trap geometries with A between 10 
and 25 and the results were unchanged. This shows that 
the thermodynamic instability of the vortex state pre- 
dicted to occur in a very elongated geometry (A > 15) 
with Cl ~ 0.75 does not prevent vortex nucleation [ p"6|| . 

In the second study that we present, the stirring fre- 
quency Q was fixed and the rotating deformation param- 
eter e was increased linearly from zero to a final value in 
the range (0.005,0.032) at a constant rate e = 0.078s -1 . 
Then, as before, we measured either a or we checked 
for vortex nucleation. In this case, since a is initially 
zero, the branch that is followed depends on which side 
of the quadrupole oscillation resonance (w/a/2) is il (see 
the dotted line in Fig. 1 corresponding to e = 0). The 
results are shown in Fig. 2, and we observe that when 
Q = 0.7 the condensate follows branch I. No vortices are 
nucleated in this case which is in agreement with the 
predicted stability of this branch |l7|. By contrast, when 
Cl = 0.75 the condensate follows the lower part of branch 
II up to a value of e c = 0.023 where branch II bends 
backwards. Beyond this point the measured value of a 
is no longer a solution of Eq. (^|) and we observe the nu- 
cleation of vortices. This situation is quite different from 
the one observed in the descending frequency ramp dis- 
cussed above. In that case, when the back bending of 
branch II is reached no vortices are nucleated and the 
condensate is observed to switch to branch I. 

In this second procedure, where Q is fixed and e in- 
creases slowly, the critical value e c at which the back 
bending of branch II occurs increases with increasing Q. 
This route can lead to vortex nucleation only if the final 
value of the ramp in e is above e c , which for a given ramp, 
puts an upper limit on the stirring frequency leading to 
vortex formation. This is well confirmed experimentally; 
more precisely, with our maximal value of e = 0.032, we 
could nucleate vortices only when the stirring frequency 
was in the range 0.71 < f2 < 0.77. 

A third route to vortex nucleation is followed when the 
stirring potential is switched on rapidly (~ 20 ms) and 
held constant for 300 ms with a fixed stirring frequency 
In particular, we found that vortices are nucleated 
along this route for ft = 0.70 and e = 0.032 (final param- 
eters in Fig. 2a) in marked contrast with the results of a 
slow ramp of e. In this case, the state of the BEC cannot 
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follow the sudden change in e. To illustrate the nonsta- 
tionary nature of the condensate we have plotted in Fig. 3 
a typical time dependence of a with this nucleation pro- 
cedure. We observe that a oscillates around 0.3 during 
the first 200 ms, and then falls dramatically to a fixed 
value below 0.1 at which point we see vortices entering 
the condensate from the border and settling into a lattice 
configuration. We should point out that the presence of 
the rotating anisotropy is not necessary beyond the first 
~ 80 ms for vortex nucleation and ordering; nonetheless, 
the number of vortices generated is an increasing function 
of the stirring time after this threshold. 
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FIG. 3. Measurement of the time dependence of a when 
the stirring anisotropy is turned on rapidly (20 ms) to a value 
of e = 0.025 and a frequency of Q = 0.7 (A = 9.2). Five images 
taken at intervals of 150 ms show the transverse profile of the 
elliptic state and reveal the nucleation and ordering of the 
resulting vortex lattice. The size of an image is 300/^m. 

In conclusion, we have studied the response of a Bose- 
Einstein condensate confined in a magnetic trap to a 
rotating perturbation created by a stirring laser beam. 
We observe an irrotational state of the condensate which 
corresponds to an elliptic cloud, stationary in the rotat- 
ing frame. We show that vortex nucleation in this sys- 
tem is related to dynamical instabilities of this irrota- 
tional state. This connection explains why the critical 
frequency originally observed in this system is notably 
larger than the one predicted by a purely thermodynamic 
analysis: when a vortex state is thermodynamically al- 
lowed but no dynamical instability is present, the time 
scale for vortex nucleation is probably too long for this 
process to take place in our system. A natural exten- 
sion of this work is the investigation of other irrotational 
steady states associated with a rotating potential of a 
higher multipole order fiofl . The exact role of tempera- 
ture in vortex nucleation remains to be elucidated. We 
believe that it has conflicting roles in the generation of an 
ordered vortex lattice (see images in Fig. 3). On the one 



hand, when the temperature is increased (while remain- 
ing below the condensation temperature), the ellipticity 
induced by the stirring is reduced and vortex nucleation 
is hindered and may even become impossible. On the 
other hand dissipation is clearly necessary to order the 
vortex lattice, and we observe qualitatively that the or- 
dering time increases with decreasing temperature. We 
plan to address these issues in future work. 
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